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1. Introduction 

Rectangular billiard perturbed by a zero-range perturbation (Seba billiard [TJ) has 
been investigated theoretically [HI3E1S1I3E1I3IH1IS1 HOI [H] and. experimentally in 
microwvave experiments with cavities [5J 1121 1101 lllj . Experimental wavefunctions 
are extracted from either the resonance shift induced by a movable perturber 
[T3l H4l [T5l [T6] or obtained by a Lorentzian fit of the reflection amplitude measured 
by a movable antenna jTTJ [TBI US] ■ However these methods work only for well isolated 
resonances. In this paper we discuss the experimental wavefunctions in the case of 
nearly degenerate states. The starting point is the general expression for the scattering 
matrix (S-matrix) for such systems, that we derived in the previous paper [10j . called 
Part I. It considered the main properties of cavities with point-like perturbations, 
following the ideas of 20 . In the derived expression the whole information on the 
antenna construction was reduced to four effective constants. We have shown that 
a lot of experimental results can be treated within the framework of the developed 
approach. We present the general expression for the S-matrix obtained in Part I in 
Section [3] of this part. 

The practical application of the derived expression for S'-matrix requires to 
compute the Green function of the unperturbed cavity as well as the so-called 
renormalized Green function [3] . In the general case the exact computation of these 
functions is by no means a trivial task. In |11) we have shown that for a rectangular 
cavity the required computation can be performed by an application of Ewald's method 
[2T1 [22] . This observation allowed us to revise the level-spacings statistics of Seba 
billiard [TJ. Ewald's method gives the representation of the (renormalized) Green 
function in terms of exponentially convergent series. This representation fits very well 
for numerical studies. 

The main purpose of the present Part II is a further application of the theory 
developed in Part I to the theoretical description of experimental data obtained in 
rectangular cavities with particular emphasis to wave functions. The experimental 
set-up and results are explained in Section [2] As it was explained above, the main 
advantage of rectangular cavities is the possibility of the exact computation of the 
S'-matrix even for considerably high frequencies. 

In Section [4] of this Part we treat our experimental findings [12] (see also figure [TJ 
obtained in the experiment with two antennas in the case of nearly degenerate states. 
We show, that while the nodal lines of eigenfunctions of a cavity are either closed or 
terminate at the boundaries, the nodal lines of the experimental patterns measured in 
this case may break inside the cavity and reduce to points in the case of degenerate 
states. This unavoidable experimental imperfection was not pointed out earlier. The 
results of Section [4] are universal and remain valid for cavities with arbitrary shapes. 

2. Experimental set-up and results 

Figure [TJb) shows a sketch of the apparatus allowing an automatic registration of 
microwave billiard wavefunctions. The billiard consists of two parts made of brass: i) 
the upper part is a plate supporting a movable antenna in the center, ii) the lower 
part is a cavity of a required shape (rectangular in our case), drilled-in into a thick 
plate. The apparatus was tested with help of a rectangular billiard with side lengths 
a = 340 mm, b — 240 mm, and a height of h = 8 mm. A fixed antenna attached to 
the lower part at the position x = 25 cm, y — 14 cm (circular white/black markers 




Figure 1. Left (a): Experimental spectra of a rectangle for different positions 
of a movable antenna. The antenna is shifted by 10 mm for each spectra in 
y direction at a fixed x = 5 cm position. Diamond/circular marker gives the 
lower/high eigenfrequency. They are obtained by fitting a double Lorentzian. 
The lines below the correspond to the modulus of the height of the resonance. 
Right: At the top (b) a sketch of the measurement set-up is shown and below 
the two resulting eigenfunctions from the double Lorentz fit using the residua are 
plotted in grey scale, (c) for the lower and (d) for the higher frequency. 



in figure [He, d)) allowed the registration of transmission spectra as well. Thus all 
components Sn, S\2, S^i, S22 of the scattering matrix could be determined by the 
vector network analyzer. Each antenna consists of a small copper wire of diameter 0.3 
mm attached to a standard SMA chassis connector and introduced through a small 
hole of diameter 1 mm into the resonator. 

Due to the presence of measuring antennas the resonances arc shifted and 
broadened. As a consequence the non-overlapping resonance approximation becomes 
insufficient in the moment where shifts and widths induced by two (movable and fixed) 
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antennas and additional broadenings due to the wall absorption are of the same order 
as the distance between the unperturbed resonances. In Figure 0Ja) the spectra for 
different y antenna positions are shown in a frequency range where the perturbation by 
the antenna is of the order of the frequency distance between neighboring resonances. 
From the residua of each double Lorentzian fit the wave functions are obtained (see 
Figure [UJc) and (d)). In these figures one sees that dark "islands" are connected 
with each other in several points via light "bridges" , which evidently break the nodal 
lines of the pattern. The appearance of these bridges and the broken nodal lines 
of the patterns are the most important experimental imperfections in the case of 
overlapping resonances. We discuss these imperfections in details in Section [5] In 
the billiard geometry described above we found several nearly degenerate states. The 
theoretical description used up to date [T2] had not been fully satisfactory, though 
describing some of the features of the obtained wavef unctions. It is the purpose of 
this paper to obtain a full understanding of the experimental pattern, and especially 
the nodal line structures. The results are also valid if the perturbing bead method is 
used to measure wave functions [13j [23l [24J [25l [26] . 



3. Sketch of the general theory 



In this Section we remind of the general formulas we obtained in Part I of the paper. 
We have found the following expression for ^-matrix of the cavity with attached 
antennas: 

-l 



A - ik 



2ikB 



A + ik A + ik 
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BC 



C 



A + ik' 



(1) 



A + ik 

where A, B, C are diagonal matrices with elements Ai, Bi, d on diagonals and index 
i denotes a number of an antenna. Here C* — Bi describe the coupling to outside and 
Ai describes the boundary conditions at the end of the i-th antenna. The matrix G 
has elements 

ft (R*;fc 2 ), 



Gij — G^ (Ri, . . . , R;,; k ) — 



G(Ri,Rj-; k 2 ), i^j. 



(2) 



Here i = \,...v are the points of antennas attachments, ^(R^fc 2 ) is the 
renormalized Green function 
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and Pi is the scattering length of the i-th antenna. Removing a global phase we can 
write the S'-matrix in the form 

-l 



5 = 1 



2ikB 
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(4) 



The resonances of the S'-matrix ((4]) are obtained from the complex zeros of the equation 



det (d- BC/(A + ikfj = 0. 



(5) 



On the theory of cavities with point-like perturbations. Part II: Rectangular cavities 5 



4. Wavefunction probing in case of overlapping resonances 



The advantage of the theory in Section |3] is its universality as compared to the n-poles 
approximation |10j . Indeed, the expression for the S- matrix (j4]) is equally valid for 
closed cavities in the cases of isolated and overlapping resonances as well as for open 
cavities. 

The S- matrix of a microwave cavity contains all available information [27j. 
Therefore any quantity of interest should be derived from the measured S'-matrix. 
Being interested in wavefunction measurements we relate below the scattering matrix 
to the experimental wavefunction. However the experimental wavefunction in its turn 
can be easily related to the eigenmode of the closed cavity only in the case of an 
isolated resonance. It was already shown [10] that for degenerate eigenvalues of the 
closed cavity one measures a sum of squares of all eigenfunctions corresponding to this 
eigenvalue. In this Section we go further and obtain a theoretical treatment of the 
experimental wavefunctions in the case of overlapping resonances. 

To probe a wavefunction one uses a movable probing antenna and measures the 
reflection coefficient. If there are together v antennas attached to a cavity, v — 1 fixed 
and one movable, the reflection coefficient at the movable antenna is given by S vv . 
Here the index v is associated with the movable antenna. In the present case there 
are just two antennas, the fixed one in the bottom plate, and the moveable one in the 
upper part. From (jj) we obtain 



Using the equalities 



\B„ 



Al + k 2 



G 



BC 



A + ik 



(6) 



det[Q] = det[K] \Q W - £ Q vi Q jv K^ 



Qzl = \ Q V u ~ J2 QviQivK^J > ( 7 ) 

where Q is an arbitrary invertible v x v- matrix, and Kij = Qij, i,j = 1, . . . , v — 1, ^ is 
invertible (y—\) x {v—1) matrix, we can rewrite the matrix element [G— BC/(A+ik)]~* 
in the following form 



G - 



BC 



A + ik 
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^(R^k^-lB^/iA^ + iky 



where 
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G^(r,R;fc 2 ) + — In 
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(9) 



is the renormalized Green function of the system with v — 1 fixed antennas [10) , 



G^- 1 Hr,R u ;k 2 ) = G(r,-R v ;k 2 )- G(R V , R t ; fc 2 )^ 1 G(R i , R,; k 2 ) (10) 

is the Green function of the system with v—1 fixed antennas and Kij = [G + D — 
BC/(A + ik)]ij, i,j = l,...,v-l. 
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4-1- Single pole approximation 

It was shown in |10) that in the single pole approximation 



G - 



BC 
A + ik 



= v- 



^(R,)/K 2 



(ii) 



where (r) is a normalized eigenfunction of a closed cavity, is the corresponding 
eigenvalue, and 



V,; = lim 



E„ - k 2 



A, + ik 



(12) 



Equality (fTTjl allows one to associate the residue of <5w with a modulo square of an 
eigenfunction of the closed cavity. The wavefunction can also be associated with a 
shift of resonance in the complex plane influenced by the motion of the v-th antenna 

u—l 

^l^i)/Vi + il>l(Ru)/Vu. (13) 



fc a l_ E \-.,.2rt>A/Tr. , ./.a/ 



The described relations between S vv and wavefunctions of the closed cavity lead to 
two experimental standard techniques to measure wavefunctions. The wavefunction 
is either related to a residue or to a shift of a resonance. However both definitions are 
based on a single pole approximation and thus possess intrinsic imperfections. Below 
we analyze both methods in the case of overlapping resonances. 



4-2. Resonance shift- 



Poles k = ky, of 



can be found from the equation 
$ -1) (R*; k) - \B V \ 2 /{A V + ik) = 0. 



(14) 



Equation (fT4)) can be easily treated if the pole structure of £g 1 (R K ; fc 2 ) is known. 
Let us find the residues of the function fg" -1 ^ (R; fc 2 ). To this end we will derive the 
convenient representation for G^^'iv, R; k 2 ). This Green function can be written in 
the form 

G {u - 1] (r, R; fc 2 ) = G (,y - 2) (r, R; k 2 ) 

G^- 2 )(r,R„_ i; k 2 )G^-V(R u ^ u -R; k 2 ) 



4^ 2) (R^ i; k 2 ) - l^.ilV^-i + ik) 



(15) 



where 

G ( "- 2) (r,R; fc 2 ) 



G(r,R;fc 2 ) - ^ G(r, R,; fc 2 )M^ 1 G(R J ■, R; fc 2 ) (16) 

is the Green function of the system with i/ — 2 antennas, My = , i, j = 1, . . . , v — 2 
(see |Appendix A[ ) . Thus the renormalized Green function (R„; fc 2 ) takes the 

form 

e^- 1 )(R,;A ; 2 ) = ^- 2) (R,;fc 2 ) 

G^ 2 ) (R,, R„_i; fc 2 )G^~ 2 ) (R„_i, R y ; fc 2 ) 



^: 2) (R,_ i; fc 2 ) - l^.il 2 /^.! + ik) 



(17) 
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Now the /z-th residue of -"-'(R^; k 2 ) can be found straightforwardly: 
Res [&- l \R„;k 2 )] k = 



G^ 2 ) (R„, R„_ i; ql)G^-V (R v _i,Rj,; g 2 ) 



(18) 



[^:; J (R,-i; P)/afc2]^ + z|S,-i| 2 /2<Zm(^-i 

where is the /K-th zero of the denominator in (|17l) . 

It can be shown that G^ -2 - 1 (r, R„_i; g 2 ) is an eigenfunction of the non- 
Hermitian operator corresponding to the cavity with v — 1 attached antennas with 
boundary condition permitting only outgoing waves. In its turnG^- 2 )(R I/ _i,r;g2) = 
(jC" - 2 )(r, R„_i;g 2 ) is an eigenfunction of the corresponding Hermitian conjugate 
operator. 



4-3. Lorentz-fit 



A Lorentz-fit of S vv obtained from the probing antenna gives residues corresponding 
V 



to the poles fc 2 of ([8]). One associates the residue 
-l 



Res 



G- 



BC 



A + ik 



1 



[^"^(R^fc 2 )/^ 2 ]* 



iM 2 



(19) 



where fc M are solutions of (|14[) . with the modulo square of the /i-th experimental 
eigenfunction. 



4-4- Small width approximation 

When Bi -> 0, i = \,...,v we can put i? 4 = in ([I3] ) -(fl9 ]) but not in fg]), since 
otherwise the resonance term in would disappear. Then all poles becomes real 
and the S-matrix elements become singular. This approximation breaks immediately 
the unitarity of the S'-matrix for a closed cavity. Nevertheless, it can be used since it 
almost preserves the residues as well as the positions of the resonances. 

Then (JTHJ) , (Tl9l) can be simplified. Indeed the Green function of a closed billiard 
with v — 1 scatterers can be written in the form 

°° q(v-2) 



G^- 1 )(r,R;fc 2 ) = X: 



»(r,R i ,_ 1 ; g g)G^- 2 )(R,_ 1 ,R;g 2 ) 



where 



N = — 



-(R*-i;d). 



dfc 2 ' ^ 

Using the equality following from (|3]) and the standard spectral sum we find 
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where tr is the area of the billiard. Thus 
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G^-^foR; fc 2 



(20) 

(21) 

(22) 

(23) 
(24) 



On the theory of cavities with point-like perturbations. Part II: Rectangular cavities 8 



From (fT5)) and (|2"5)) we conclude that the renormalized Green function l > takes 
the form [TT] 



(R; fc 2 ) = CR5 0) + E /2(R) {-l\ 2 4) , (25) 

where k ~ are the solutions of the equation 

^: 2) (R,_ i; fc 2 )=0 (26) 

and 

f,{v) = G^{v,K u - l -ql)/N, (27) 

are normalized eigenfunctions of the system with v — 2 scatterers. 

From (IT9l) and (|24| we find that the pattern associated with a residue is 
proportional to 

ff) (28) 

/ a d»r [G(-D(r,R;fc2)] 2 

where k = k^ are the solutions of the equation 

47 1) (R,;fc 2 )=0. (29) 

5. Experimental patterns in the case of nearly degenerate states 



The universal form of the representation (|25p allows one to analyze main features of 
experimental wave-patterns equally for any number of introduced antennas. Especially 
simple analysis can be performed if the experimental patterns are associated with 
changes of resonance positions while the probing antenna is moving inside a cavity. 
We perform this analysis for the case of nearly degenerate states. 

In |10] we have shown that in the case of a degenerate state the observed pattern 
is a sum of squares of all wavefunctions forming the orthogonal basis in the subspace 
corresponding to the degenerate eigenvalue. Thus an experimentally measured pattern 
has no longer a simple relation to an individual wavefunction of the unperturbed 
system. 

What is happening in the situation with nearly degenerate states? It is clear that 
when a shift of a resonance due to a presence of a measuring antenna is comparable 
with the level-spacings, one can not treat an observable wave-pattern as a square 
of an isolated wavefunction. One can expect at least that the nodal lines of both 
eigenfunctions remain in the experimental patterns, since an eigenstate can not be 
excited by an antenna moving along its nodal line. However this expectation fails. We 
show below that in this situation the measured wavefunctions are very peculiar: they 
cannot only be related to isolated wavefunctions, but the nodal lines of these patterns 
are broken. 

For our purpose it is sufficient to consider the renormalized Green function (|25|) 
in two poles approximation. Substituting it into (|14p with B v = we obtain 

|S> + |5| + D(R)= 0, (30) 
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where R = is the position of the probing antenna, z^ = g 2 , 

/i( R) = G^PM^-i;^) , (31) 



D(R) 



'^(Rjfe 2 ) /l(R) /I(R) 

^" ' zi — k 2 Z2 — k? 



(32) 

fc 2 = (zi+z 2 )/2 

We assume that |D(R)| is large enough to neglect the influence of all other poles. 
Under this assumption (|30| has two positive solutions k 2 = k\,k\. 

Consider first the possibility of the degeneracy of two solutions k\ — z\, k\ = Z2 of 
(|30|) . The condition k\ (R) = k 2 (R) seems to define a curve in the plane of the cavity. 
However from the form of (I3TJ1) one can see that there are always two different solutions 
except for /i(R) = and fc 2 (R) = z\ (analogously if /2(R) = and fc 2 (R) = 22)- 
Thus the degeneracy can occur only if a scatterer is placed somewhere at a nodal line. 

Since the degeneracy of eigenvalues can occur only at isolated points, in the 
vicinity of these points eigenvalues should split. Consider the nodal line r = R(t) of 
the function / x (r): /i(R(r)) = 0. Along this line j3flj reads / 2 (R(t))/[z 2 - fc 2 (r)] + 
£)[R(t)] = 0. Let us assume that k 2 (r) — z\ for some value of the parameter r = tq. 
At the point Ro = R(tq) the eigenvalue z\ is twofold degenerate. 

Let us now consider the line r = R(t) +£(t), where |£(r) > is small enough. 
The solutions k\ 2 of ([30|) along this line read 



,2 _ £>(R + CKfi + z 2 ) + A 2 (R + g) + /|(R + j)± VP(R + Q 

kl ' 2 ~ 2£>(R + £) ' [66} 

£>(R) = [d(R)( Zi - z 2 ) + / 2 (R) - /|(R)] 2 + 4/ 2 (R)/|(R) (34) 

where plus and minus signs correspond to different branches, since T> > 0. Far from 
the point R = Ro we can write 

,2 „ gCgOCgi + * 2 ) + /l(R) ± I J>(R)(*i - gg) - /|(R)| 

fcl < 2 " 2~d(r) (d5j 

or 

,2 ^ ; 2 , /f(R-) -r /K-^O . n/ -Tj\ ^ n 

k x ~zi, fc 2 = z 2 + , if h-D(R)>0, 

L>(R) z 2 - zi 

+ fc 2 2 = ^i, if ^ R 1 + ^(R)<0. (36) 

D[R) z 2 - zi 

In the vicinity of the point R = Ro the expression 

+ D(R) (37) 

Z 2 — ^1 

changes its sign. Thus in the vicinity of this point along the line r = R(r) + £(t) an 
avoided crossing takes place. 

Figure [2] illustrates what may happen in a wavefunction measurement in the 
case of nearly degenerated eigenfunctions. The calculation has been performed for 
a rectangle with a side ratio of 1/1.01, i.e. very close to the square. Figures (a) 
show the eigenfunction with quantum numbers (2,3) and (3,2), respectively, for the 
rectangle without perturbation by an antenna. Due to the small side mismatch there 
is an energy splitting of the the eigenstates of the rectangle. 

The avoided crossing allows us to separate properly the patterns corresponding 
to lower and higher frequencies (see figure [2] (c)). It also shows that a nodal line of an 
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Figure 2. The expected experimental patterns (squares of experimental 
wavefunctions) in the case of two nearly degenerate states of a rectangular 
billiard. The side ratio d x /d y of the rectangle is 1/1.01. Row (a) shows a 
pair of eigenfunctions of the rectangle with quantum numbers (2,3) and (3,2), 
respectively. The next two rows show the evolution of the experimental patterns 
corresponding to the higher and to the lower eigenvalue, respectively, for two 
different position independent perturbations D=l (b) and 0.1 (c). The nodal 
lines are inserted into the pattern as black solid lines for better visuality. In the 
left column all nodal lines remain unperturbed, in the right one they become 
interrupted (b,c) and tend to points for (d). The right column show shifts of 
two resonances (rr-axes) along the left vertical nodal line in Figures (b, c) (solid 
curves) and a line close and parallel to it (interrupted curves). In figure (d) the 
sum of figures (a) is plotted. 
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unperturbed eigenfunction may become interrupted (figure[3](b,c) center column), and 
the "bridges" connecting nodal domains through the nodal line may appear. These 
bridges are expected to appear in the case when a shift of a resonance due to an 
antenna or a perturber is larger than the level-spacing. 

On the contrary, the structure of nodal lines in the first column in figure [2] remains. 
This is clear since for both resonances the shifts due to the presence of a scatterer 
has always the same sign. When the influence of the perturbation increases (|D(R)| 
decreases), the pattern keeping the nodal line structure of the original unperturbed 
state is concentrating in vicinities of nodal lines of the the other unperturbed state. 
The other pattern tends to the sum of squares of two unperturbed wavcfunctions (see 
figure [2](d)), i.e. tends to the pattern expected for a degenerate state [TO], 

One can see that measured nodal lines for the case of a single scatterer are 
continuous if the coupling of a scatterer is weak as compared to level spacing. 
For a stronger coupling nodal lines become interrupted and turn into the points of 
intersection of nodal lines of f\ (r) and fa (r) . 

5.1. Experimental verification 

Now we have all the insights to treat the experimental findings presented in Section [5] 
The experiment was performed on the rectangular billiard with sides d x = 34cm and 
d y = 24cm. The square of the side ratio (d y /d x ) 2 = (34/24) 2 ~ 2. The eigenvalues of 
the unperturbed billiard are 



It is easy to show that n 2 + 2m 2 can be degenerate. An example is n\ = 5, mi = 
1, n 2 — 3, m 2 = 3. 

Figure [3] illustrates the comparison of numerical simulations with experimentally 
measured patterns corresponding to the eigenstates with quantum numbers n = 
5, m = 1 and n = 3, m = 3. In numerical simulations we assumed for simplicity 
that D±, Z?2 do not depend on R. We also assumed that the coupling strength 
of both antennas is negligible, thus the resonances can be found from zero-width 
approximation . 

In the figures EK),|3b) one sees the original eigenfunctions of the cavity. In figures 
Et),[3jl) we neglected the splitting of eigenvalues as well as the influence of the movable 
antenna. Figures Ek), |3f) illustrate the numerical simulation with fixed and movable 
antennas. Plotted values are the solutions k 2 of (|T4"|). Finally figures [5g),[3h) represent 
the experimental data. One can see that there is a reasonable coincidence between 
experimental figures and the numerics. 

This work was supported by the Deutsche Forschungsgemeinschaft via an 
individual grant and the Forschergruppe 760: Scattering systems with complex 
dynamics. 

Appendix A. Representations of the Krein's resolvent 

Krein's resolvent G^\r,'R; k 2 ) of the billiard with v antennas reads 




(38) 




(A.1) 
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Figure 3. Calculated distribution \ip nm \ 2 . White and black points show the 
position of the fixed antenna x = 25cm, y = 14cm. a,b) Without antennas: 
IV'iml 2 : a) n = 5, m = 1, uigi = 2. 29113GHz; b) n = 3, m = 3, = 
2.29348GHz. Here and further fc^m" 1 ] ~ 0.2096w[GHz] . c) |G(r, R, fc 2 )| 2 , 
u) = 2. 28706GHz. d) The combination of eigcnfunctions a) and b) equal to 
zero at the antenna point. e,f) Simulation of the experimental wavcfunctions 
in the billiard with two coupled antennas using the Ewald's method 11, and the 
standard procedure of the Lorenzian fit 1 171 . We used the following parameters: 
A 1 = A 2 = 0.3dm- 1 , B 1 = B 2 = 2.5dm- 1 / 2 , D x = -25, D 2 = -65. g,h) 
Experimentally measured wavcfunctions: g) u = 2.298GHz, h) lu = 2.307GHz. 
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From the physical point of view it is obvious that instead of the unperturbed billiard 
corresponding to the free Green function G(r, R; fc 2 ) one can consider the billiard 
with m < v antennas corresponding to the Green function G t - m ' (r,H; k 2 ) as the 
unperturbed one and the rest of v — m antennas as the perturbation. Then (|A.1[) 
can be rewritten in the form 



G M (r, R; k 2 ) — G (m) (r, R; k 2 ) - ^ G (m) (r, R i; k 2 ) 

»J=1 



Jm) 

Hi 



G (m) (Rj,R; fc 2 ),(A.2) 



where Q 



(m) 



[G^+D-BC/(A + ik)] l3 , i,j 1 v, 



G 



(m) 



Gy"^(Ri, . . . , R„; k 2 



4 m) (R i; P), i = j, 
G^HR^R^k 2 ), i^j. 



47 ) (R;fc a )= lim 

r > R 



G (m) (r,R; k 2 ) + — In 
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Ir-R| 



(A.3) 
(A.4) 



The formal proof follows from the equality 

Qn Q12 \ 
Q21 Q22 / 

= ( Qn + QuQi3 [Q (m) ] ~ l QaiQn 1 -Qn Qia [Q M ] ~* ^ 

valid for any invertible matrices Qn, Q22- Here Qn is m x m, Q22 is (i/ — m) x (v — m) 
matrix, Q< m ) = Q 22 - Q2iQnQi2- Substituting (|A.5|) into (|A.1|) we find 

G (l/) (r,R; fc 2 ) = G(r, R; fc 2 ) - (r,m|Qi 1 1 |R,m> 



(A.5) 



-((*,*> 



(r, m|Q 11 1 Qi2 



Q 



(m) 



(|R,i/-m) -Q 2 iQn\R,m}\ (A.6) 



where (t.v — 



G(r,R i; fc 2 ), (r,m\j = G(r,R f ,k 2 ), \R,v - m) % = G(R i; R; fc 2 ), 
|R, m)j = G(Rj, R; fc 2 ), i = 1, . . . , z/ — m, j = 1, . . . ,m. Taking into account equalities 

'(r,i/-m| - (r,m|Q n 1 Qi 2 ) i = G^ (r, R i5 fc 2 ), 

|R, ^-rn) - Q 2 iQ n 1 |R, m)^ = G (m) (Rj , R; fc 2 ), (A.7) 

we conclude that (|A.6[) is identical with (|A.2[) . 
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